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Abstract 

The generalized Stieltjes transform (GST) is an integral transform that depends 
on a parameter p > 0. In previous work a convenient form of the inverse trans- 
formation was derived for the case p = 3/2. This paper generalizes that result to 
all p > 0. It is a well-known fact that the GST can be formulated as an iterated 
Laplace transform, and that therefore its inverse can be expressed as an iterated 
inverse Laplace transform. The form of the inverse transform derived here is a 
one-dimensional integral that is considerably simpler. 



1 Introduction 



In [T] we encountered the integral equation 

POO 

9(0= (C 2 + ^T 3/2 /(aO 4< |argC|<vr/2, (1) 
Jo 

where g(() was a known function, and we needed to solve for /(/x). That paper proved that 
if g has suitable analyticity and asymptotic properties, which were satisfied for the specific 
function of interest, then 

f(ji) = -V- [ cos6 g(-ifMcos6)dO. (2) 



7T 

This result is reviewed in 

By a change of variables, eq. (J2J) can be rewritten in the form 

/(A0 = £ 7= g(0C (3) 

Then letting y = fj 2 and z = £ 2 , as well as F(y) = ^-f(fi) and G(z) = g((), we obtain 

/>oo 

G(z) = (y + z)-^ 2 F{y) dy, | argz| < vr (4) 



and 

W 4tTZ 4 ^T~z 1 ; 

The contour C y starts and ends at — y and encloses the origin in the counterclockwise sense. 
It is convenient to have a y-independent contour, so letting z = wy we obtain 

F to )= 1 VS/^L*., (6) 

where C is a contour that starts and ends at the point w = —1 enclosing the origin in the 
counterclockwise sense. It could be chosen to be the unit circle, for instance. Note that it 
is not really a closed contour, since G(z) has a branch cut running along the negative real 
axis. An integration by parts allows this to be rewritten in the form 

F(y) = -^-y 3/2 [ VTT^G'(yw)dw. (7) 
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In this paper, we will formulate and prove a one-parameter generalization of the preceding 
result. Specifically, we claim that the integral transform 

POO 

G(z)= {y + z)-PF{y)dy, |argz|<7r (8) 
Jo 

1 



has as its inverse transform 



f(v) = ~ v p /(i + wy-'G'iyw) dw 



(9) 



Jc 



This integral converges for p > 0. The special case discussed above corresponds to p = 3/2. 
The counterpart of eq. (JBJ), obtained by integration by parts, namely 



is less general, since it is only valid for p > 1. This is allowed in the p = 3/2 case, of course. 

We should note for the record what is being assumed about F and G. Specifically, F(y), 
which is only defined on the positive real axis, is allowed to be an arbitrary distribution (or 
"generalized function"). We require that there exists a number a with < a < p such that 
| Jy* y a ~ p F(y) dy\ is bounded by a number independent of y\ and y 2 for all < y\ < y<i- The 
function G(z) is then holomorphic throughout the cut plane | arg,z| < 7r, and there exists a 
positive real number (3 such that |2r G(2)| is bounded at infinity. 

The p = 1 case of eq. (|HJ), known as the Stieltjes transform, is discussed in Widder's 
classic treatise on the Laplace transform The only other cases considered by Widder are 
positive integer values of p, which are related to the p = 1 case by differentiation. Following 
[3], we refer to the case of arbitrary p as the generalized Stieltjes transform (GST). (In ref. [H] 
it is called a Stieltjes transform of index p.) The formula for the inverse GST in eq. does 
not seem to have been found previously. Certainly, it does not appear in [3], 0, or [Bj. 

Much of the literature on the GST is concerned with the asymptotic behavior of G(z) for 
large \z\ We will not address that topic here. As it happens, in pQ we were interested in 
deducing the asymptotic behavior of F associated with a given G. 

In the special case p = 1 the transform in eq. (JSJ) reduces to the Stieltjes transform 



Since the Stieltjes transform is well understood, this case provides an instructive test of 
the proposed inverse transform. Setting p = 1 in eq. gives an expression that can be 
integrated explicitly to give 



e^o+ 2iri V / 

Given the stated analytic and asymptotic properties of G(z), it is a simple consequence of 
Cauchy's theorem that this is the correct solution of eq. (fTTj) for y > 0. 
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(10) 




(11) 




(12) 



In order to convince oneself that eqs. (jSJ) and (JSJ) are plausible for all p > 0, it is instructive 
to consider a simple example. Specifically, if one chooses 

F(y) = y v ~\ (13) 

then the integral in eq. (jHI) converges for < v < p and gives 

G(z)=2T-"B(u,p-u), (14) 

where B(u,v) is the Euler beta function. It is straightforward to verify that this pair of 
functions also satisfies eq. (JHJ). 

2 Other versions of the inverse transform 

The change of variables z = wy allows us to rewrite eq. © in the alternative form 

F(y) = -^- [ (y + zy- 1 G'(z)dz,. (15) 

2™ JCy 

As before, the contour C y starts and ends at the point z = —y, encircling the origin in the 
counterclockwise sense. If p > 1, an integration by parts brings eq. ()15j) to the form 

F(y) = ^-(p-l) [ (y + zy- 2 G(z)dz. (16) 

2711 JCy 

Since this is only well-defined for p > 1, eqs. Q and ([To]) are more general than eq. (fTB^ . 
Let us define the quantity that appears on the right-hand side of eq. (fT2|) to be 

A(t) = — lim (G(-t - ie) - G(-t + ie)), t > 0. (17) 

We showed that F(y) = A(y) when p — 1, but this is not the case for other values of p. By 
shrinking the contour C y down to the cut, eq. (|T6*j) takes the form 

F(y) = (p-1) f\y-ty- 2 A(t)dt. (18) 

In similar fashion, eq. (|15|1 gives rise to 

F(y)= !\y-ty- l A\t)dt. (19) 



o 



However, these formulas are only correct if the behavior of G(t) near the origin is such that 
these integrals exist. The contour integral versions of these formulas are more general, since 
they do not have this restriction. 
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Equations (|18|) and (J19)) have the structure of Abel transforms. The inverse Abel trans- 
form is well-known and can be used to give a formula for the discontinuity across the cut, 
A(t), in terms of the original (generalized) function F(y). A version that is suitable if p < 2 
and F(0) = is 

A(t) = !\t - yf-PFXy) dy. (20) 

71 (1 - P) Jo 

This gives the known result for p — 1, namely A(t) = F(t). It can be checked for the simple 
example F{y) = y v ~ x discussed earlier. 



3 Proof of the main result 

Let us review how the GST is related to the Laplace transform, for which we use the following 
notation 

/■oo 

C X [F]= I e- xy F{y)dy. (21) 

Inserting the identity 



into eq. (jSJ) gives 



o 



( y + z )-p= x^e-^+^dx (22) 

r(p) Jo 



1 [CO 1 

G ( z) = T\pjJo xP ~ le ~ xzc ^ dx = Y^ c *[ xP ~ lc ^~\- ( 23 ) 

This recasts eq. (jHJ) as a Laplace transform of a Laplace transform. In particular, setting 
p = 1, this gives the well-known result that the Stieltjes transform is the square (in the 
operator sense) of the Laplace transform, i.e. S = C 2 . 
Equation (|2^|) implies that 

C x [F]=T( P )x 1 ^C- 1 [G}, (24) 

where £ _1 denotes an inverse Laplace transform. A second inverse Laplace transform gives 
the formal inversion of eq. ()23|) 

F(y) = T(p)C y 1 [x 1 -»C- 1 [G]\. (25) 

This can be made very explicit by using the standard contour integral realization of the 
inverse Laplace transformation [Sj. The new claim is that eq. (j25j) can be simplified to take 
the form of eq. (JHJ). 
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Let us now carry out some similar manipulations of eq. (JOJ). Substituting 

POO 

G'(yw) = - dte- tyw tC- l [G] (26) 
Jo 

into eq. (JSJ) and taking the Laplace transform of both sides recasts eq. Q in the form 

C X [F] = T(p + 1)— / dw(l + w) p - 1 [ dt(x + tw)- p -HC^ l [G\. (27) 

27U Jc J 

Comparing eqs. (f2~2j) and (|2*7j) . we see that their equivalence requires the identity 

[ dw(l + wy-\x + tw)-"- 1 = x l - p 8(x - t), (28) 
2m Jc 



or equivalently 

P 

2m 



J dw{\ + w) p -\x + tw)-?- 1 = x- p 5(x - t). 



(29) 



Another way of understanding the necessity of eq. (}2*9*j) is to consider the special case 
F(y) = 8(y — t). In this case the GST is G(z) = (t + z)~ p . The inverse transform eq. Q 
for this choice of G(z) corresponds precisely to eq. (|2*H|) . This remarkable equation, which is 
required to hold for x, t > 0, is the heart of the matter. Its proof, which is rather nontrivial, 
is presented in the next section. 

4 Representation of a delta function 

This section proves the key formula, namely eq. ()29|). Consider the left side of eq. (J2*§|) 

x (x,t) = ^ [ dwil + wy-'ix + tw)-?- 1 , (30) 
2m J c 

where x,t > and C is the contour shown in fig. 1. By making the change of variables 
w — > 1/w, it is easy to prove that Xp( x ,t) = Xp(t, x )- F° r t > x the singularity structure is 
shown in fig. 1(a), and the contour can be pushed off to infinity, giving zero for the integral. 
For t < x, the contour encloses no singularity (see fig. 1(b)), so the result is again zero, as 
required by x <-> t symmetry. 

Since Xp{ x ^) vanishes for t ^ x, it must be some sort of distribution concentrated at 
t = x. If we assume that it is proportional to a delta function, i.e. f(x)6(x —t), it is easy to 
derive f(x). Integrating over t from to y, gives 

f(x)9(y -x) = / —(1 + wY- 1 x- p -(x + yw)- p . (31) 

2m Jc w L J 
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(a) (b) (c) 

Figure 1: The contour C starts and ends at the point w = — 1 enclosing the origin in the 
counterclockwise sense. Three possible branch cut configurations are depicted. 

The first term on the right is an elementary contour integral and gives x~ p . The second term 
is also easily evaluated by Cauchy's theorem. It gives — x~ p if x > y and zero if y > x. Thus 
we deduce that f(x) = x~ p , and hence 

X P {x, t) = x- p S(x -t)= r p 5(x - t). (32) 

Eq. (|3*2*|) is the correct result, but the derivation given above is not rigorous. It assumes 
that the distribution is proportional to a delta function, and that it does not involve any 
derivatives of delta functions. A more careful analysis that is sensitive to such terms if they 
are present involves checking the proposed answer by integrating both sides against the test 
function e~ ty , i.e., comparing the Laplace transform of both sides of the equation. Thus we 
need to show that 

/■oo 

dt X p(x,t)e- ty = x- p e~ xy . (33) 



o 



Evaluating the Laplace transform of eq. (|30p gives (after some simple manipulations) 



o 27 ™ jc 



dt X p(x, t)e' ty = ^—y p dw(l + w) p ~ l I{xyw), (34) 



where 



I( u ) = e u e- v v- p - l dv. (35) 

J u 

In order to reexpress this function in a more convenient form, we first note that 

I'(u) = I(u) - u- p -\ (36) 

Let us now assume that p is not an integer, which is the case of most interest. This 
equation is then solved by an expression of the form 

I(u) = h(u)+u- p I 2 (u), (37) 
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where h(u) and I^u) are regular at u = 0. The first term is h(u) = T(—p)e u . Its 
contribution to eq. (J3~4^ is zero, since there is no singularity inside the contour. The function 
that matters is I2, which satisfies the differential equation 

ul' 2 {u) -(u + p)I 2 (u) + 1 = 0. (38) 

Substituting a power series expansion, 

00 

I 2 (u) = Y,CnU n , (39) 

n=0 

one obtains the recursion relation (n + 1 — p)c n+ i = c n . Thus, since c = -^(0) = 1/p, we 
conclude that 

C = j (l - P) y (40) 

pT(n + 1 — p) 

To determine the contribution of each term in eq. (|3~9"j) to eq. (jHH), we need to evaluate 



/ dw(l + w) p - l w n ~ p = i —B(n -p+l,p) = '—. (41) 

2m J c it ' pc n n\ 

Combining these results, the c n factors cancel, and we learn that 

{-xy) n 



poo . . 

/ dt Xp (x,t)e~ t y = x- p J2 

JO „_ n 



n=0 



x- p e~ x \ (42) 



which is the result we set out to prove. Even though this derivation needs to be modified 
when p is a positive integer, the result clearly is valid in that case as well. This completes 
the proof of eq. (}3*3*j) and hence of the inverse transform eq. (|9*jl. 

5 An alternative derivation 

In this section we present a simpler, though less general, derivation of the inverse transform. 
It follows from Cauchy's theorem and the required analytic and asymptotic properties of 
G(z) that 



JO Z + 1 

where A(t) is defined in eq. ()17|) . The validity of this formula requires that A(t) is not 
too singular as t — > 0, so that the integral exists. This is a significant restriction, since if 
F ~ y u ^ 1 with v > for small y, then G(z) ~ z v ~ p for small z. Thus one would need that 
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p < v + 1. No such assumption has been made previously, which is why the derivation in 
this section is less general. 

Equation (f£!j) is precisely a Stieltjes transform, G = S[A]. We noted earlier that this is 
an iterated Laplace transform, S = C 2 . Therefore, 

C-'IG] = C X [A}. (44) 

Comparing this with eq. (|24|h which we obtained from eq. (jHJ), we learn that 

C x [F} = T(p)x 1 - p C x [A]. (45) 

Now using CJjf" 1 ] = T(v)x~ u , this becomes 

C x [F} = (p-l)C x [tP- 2 ]C x [A}. (46) 

By the convolution theorem, this implies that 

F{y) = (p-1) f\y-ty- 2 A(t)dt, (47) 
Jo 

which is eq. fllSj) . The alternative form in eq. extends the range of validity for p, but is 
even more restricted in its requirements for the behavior of A(t) at the origin. As we pointed 
out in sect. 2, there is no such issue for the corresponding contour integrals. 
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Appendix: Examples 



In this appendix we list pairs of functions that are related by the generalized Stieltjes trans- 
form and its inverse given in eqs. (JHJ) and Q. The two simple examples in the table were 
discussed in the introduction and in sect. 2. Four pages of additional examples can be found 
in Sect. 14.4 of [3]. One of them is given on the last line of the table. The example in PQ, 
which is specific to p = 3/2 (though it can probably be generalized), is not contained in 



F(y) = -^-yfJ c (l + w)P- 1 G'(yw)dw, p >0 


G(z) = J °°(y + z)-PF(y) dy 


y v ~ x , < v < p 


B(v,p-v) z u -f 


5(y-t), t>0 


(t + z)-p 


y v ~ l {l + y)- x , 0<u<p + X 


B(v, p + X-u) z v -p 2 F 1 (u, \;p + \;l-z) 



Table 1: Pairs of functions that are related by the generalized Stieltjes transform. 
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